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; ^ La presente invention concerne uiie' metho.de pour, former, .a partir de' 
/ ^ obtenues par exploration d'une zone d'un rnilieu heterogene, un modele representatif 'de la \ 

distribution dans la zone d'une grandeur physique affranchi, (au moins partiellement) de.la [[ 
10 presence de bruits correles pouvant entacher les don^^ . \ , • . \ -, , 

La methode s'applique par exem^ 
dans une zone souterraine. -, ' ■ , / ■ ' " ;. ; i - . y] _ 

ETATDELA TECHNIC l' • 

: La deinarche- consistant des . mesur^^^ 

15 experimenta^les, s'est developpSe dans quasinieht -tous le?, (loni'aih^^ la 
technologie. Urie telle; approche est connue sous diffefentes appellations : methode desv 
moindres caiTes pour V estimation de parametres, pour la resolution du resolution/ du ' 
probleme inverse. Pour une bonne presentation de cette approche dans le contexte des 
geosciences, on pent serefererpar exemple a: 

20 - Tarahtola, ' A^.-i « Inverse Problem Theory :' Method for. Data Fitting' and Model 
Parameter Estirnation >>. Ehevier, Amsterdam, 19 87 ;- 

Notons que le terme « moindres carres » fait reference au can-e de la norme dans 
Fespace des donnees pour quantifier recart entre la reponse d'un modele (laquelle est 
r image du modele par un operateur de modelisation prealablement choisi) et les donnees, 
25 fonction cout qu'il s'agit de mmimiser pour resoudre le probleme. L'utilisatjon du caire de- 
la norme pour definir la fonction cout ne constitue qu'une commodite pratique mais n'a 
fondamentalement rien d'essentiel. D'ailleurs de nonibreux auteurs utilisent, pour des 


raisons diverses, une autre definition de la fonction cout mais cette definition reste basee 
sur r exploitation de la nomie, ou d'une semi-norme, dans I'espace des donnees. Enfin 
nous disposons d'une grande latitude pour choisir la norme (ou la semi-noiTne) dans 
I'espace des donnees (nous ne sommes en aucune fagon contraints a utiliser la norme 
5 euclidienne). Dans le cas de donnees bruitees la solution peut dependre sensiblement du 
choix fait a ce stade. Pour de plus amples developpements au sujet de ce probleme, on peut 
se referer par example aux publications suiv antes : 

- Tarantola, A. : « Inverse Problem Theory : Method for Data Fitting and Model 
Parameter Estimation ». Elsevier, Amsterdam, 1987 ; Renard et Lailly, 2001 ; Scales et 
10 Gersztenkom, 1988 ; Al-Chalabu, 1992. 

Les resultats de mesures sont souvent entaches d'en-eurs. A ces erreurs de mesures 
se rajoutent, quand on veut confronter T expert mentation avec des resultats de 
modelisation, des bruits de modelisation : les niodelisations ne sont jamais parfaites et 
correspondent done toujours a une vue simplifiee de la realite. Nous ferons done apparaitre 
15 dans la suite le bruit comme la superposition : 

de composantes non correlees (bruit blanc par exemple) ; 

de composantes correlees, ce qui veut dire que Texistence de bruit sur. un echantillon 
de mesure se traduira par Texistence d'un biiiit de meme nature sur certains des points 
de 'mesure voisins ; typiquement les bruits de modelisation font partie de cette 
20 ^ categoric. 

Lorsque des bruits correles entachent les donnees, la qualite du modele estime par 
. resolution du probleme inverse peut s'en trouver gravement affectee. Comme nous I'avons 
dit aucun operateur de modelisation n'est paifait. C'est.donc toute la communaute des 
personnes impliquees dans Tidentification de parametres decrivant un modele qui voit sa 

25 . demarche entravee par I'existence de bruits correles. Parmi ces personnes, les praticiens de 
Texploration sismique figurent parmi les plus concernes : leurs donnees ont en effet un 
rapport signal sur bruit mediocre voire tres mauvais. Voila pourquoi les techniques de 
filtrage des .bruits correles constituent une paitie impottante des logiciels de ti'aitement de 
donnees. sismiques. Les techniques les plus classiques explqitent une transformation 

30 (transformation de Fourier par exemple) ou signal et bruit se trouvent localises dans des 
regions de. I'espace differentes permettant ainsi une separation du signal et du bruit. Pour 
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un presentation generale des methodes classiques d'elimination du bruit sur des donnees 
sismiques, on peut se referer au livre de Yilmaz : 

- Yilmaz, O. 1987: "Seismic data processing", Investigation in geophysics Society 
of Exploration Geophysicists, Tulsa, 1987. 

5 Ces techniques ne sont cependant pas parfaites : elles supposent que Ton a trouve 

une transfonnation permettant la separation complete du signal et du bruit. Dans ce 
coritexte, les bruits correles apparaissent comme particulierement genants car ils peuvent 
etre difficiles .a separef .du signal (qui, liii aussi,, est correle) et , peuvent etre de fortes 
. . amplitudes. On est ainsi souv.ent conduit a la difficile gesti on du compromis suiyant : ou 

10 bien on preserve le signal maiis il reste un fort residu de bruit,^ on bien on elimine le.'bruit 
mais oh a, du coup, distofdu ie signal. Ces techniques de filtrage peuvent etre mises en 
oeuvre prealablement a la resolution du probleme inverse : elles constituent alors un 
pretraitement applique au donnees. La qualite de la solution du probleme inverse depend, 
alors fortement de raptitude des filtres a eliniiner le bruit ^ans distordre le signal. 

15 L-approche,introduite par Ne^^ \ 

, - -Nemeth T, et '4 (1999^^^ Migration- of Incomplete reflection . Data"; ; 

;-v / -Geophysics, 64, 508-221, \: ^. _ ■ , 

^ constitue une avancee importante pour rinyefsibii die donnees perturbees. par des bruits 
correles de forte amplitude. Ces auteurs proposent d'eliminer le bruit par resolution d'un 
,20 probleme inverse : ,ayant defini I'espace des bruits correles comme I'espace image d'un 
' ^'v • espace vectoriel B (espaee de fonctions generatrices de bruit) par un operateur lineaire T et 
■ ,]e signal comme Tespace image d'un espace vectonel M (espace des modeles) par un 
■ ' • operateur de modelisation F (qu'ils supposent lineaire), ils lecherchent le signal dans F(M) 
.-et le bruit. dans T(B) dont la somme est la plus proche possible (au sens de la noime 
25- euclidiehne ou, en continu, de la norme L^) de la donnee . rhesur6e. Cette technique 
constitue une avanc6e irhportante dans la mesure pu 'elle a permis 1 'elimination (ou du 
moins une reduction tres substantielle) de bruits correles de forte amplitude (ondes de 
surface pai' exemple) difficileinent separables du signal par les techniques conventionnelles 
de fikrages, Cependant, de Taveu merae des auteurs, ces performances se payent par une 
30 ^augmentation d'un ordre de grandeur du temps de calcul requis pour la resolution du 
probleme inverse classique c'est a dire sans rechercher la composante correlee du bruit. Par 


. ailleurs le resultat foumi par la methode s'avere extremement sensible a toute inexactitude 
introduite lors de la definition de I'operateur T : c'est la la contrepartie ineluctable de la 
grande aptitude de la methode a discriniiner signal et bruit. 

LA METHODE SELON L'INVENTION 

5 La nnethode selon I'mvention permet d'estimer, a partir de donnees obtenues par 

exploration d'une zone souterraine, un modele representatif de la distribution, dans la zone, 
d'au moins une grandeur physique, ce modele etant largement affranchi de la presence de 
bruits correles pouvant entacher les donnees. La methode comporte essentiellement les 
etapes suivantes : 

10 a) acquisition de mesures renseignant sur certaines caracteristiques physiques de la 

zone en suivant un protocole experimental defini ; 

b) specification d'un operateur de modelisation, qui, a un modele de chaque 
grandeur physique associe des donnees synthetiques constituant la reponse du modele, les 
mesures et les donnees synthetiques appartenant a un espace de donnees ; 

15 c) pour chaque bruit correle repere par un indice j allant de 1 a J, choix d'un 

operateur de modelisation qui associe un bruit correle a une fonction generatrice du bruit 
appartenant k un espace prealablement defini de fonctions generatrices du bruit ; 

d) specification d'une norme ou d'une semi-norme dans I'espace des donnees ; 

e) specification, d'une semi-norme dans Tespace de fonctions generatrices de bruit 
20 pour laquelle chaque operateur de modehsation du bruit etablit sensiblement une relation 

isometrique entre Tespace des fonctions generatrices du bruit et I'espace des donnees; 

f) definition d'une fonction cout quantifiant Tecart entre les mesures d'une part et la 
superposition de la reponse du modele et des bruits correles associes aux fonctions 
generatrices de bruit d' autre part ; et ' 

25 g) ajustement du modele et des fonctions generatrices du bruit en minimisant la 

fonction cout, par une methode algorithmique tirant partie, des proprietes d'isometrie des 
■ operateurs de modelisation du bruit. - " . 
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Suivant un premier mode de mise en oeuvre, on recherche la distribution en fonction 
de la profondeur de rimpedance acoustique dans le milieu, les bruits correles affectant les 
donnees sent des ondes de tube identifiees chacune par des paiametres caracterisant leur 
propiagation, les donnees mesur^es sont des donnees de type PSV en ulilisant des capteurs 

5 adaptes a detecter le deplacement des particules dans le milieu en reponse a une excitation 
sismique localisee, remplacement des capteurs, la duree d'enregistrement et les points 
d'echantillonnage en temps etant definis, et I'operateur de modelisation choisi associe les 
; donnees synthetiques a une distribution d'impedance acoustique en .fonction de la., 
profondeur evaluee en temps, de parcours et a la contrainte verticale miesuree en fonction 

10 du temps a la profondeur du preihier capteun . . / , . 

Comme fonction cout quantifiant I'ecart entre les mesures d'une part et la' 
superposition de la repofise du modele et des braits correles associes aux fonctidns 
generatrices de bruit d'autre part, on choisit par exemple le carre de 4a-semirnorme de cet.. 
ecart dans Tespace des donnees; . ' , i '\' 

15 / Suivant. un exemple de realisation, rajustement* du modele let ' des fonctions',,.; 

> generatrices du bruit est. pbtenu par exemple par ■uiie methods de relaxation pour eliihiner- 
les inconnues correspondant a. chaque fonction gerieratrice de bruits cprteles, cette methode; ' 
de relaxation etaht niise e^ geiivre a Tinterieur des iterations = d'un;a^^^ quasi- ;^ 

Nemon pour le cklcul du mod \ . ; ; " ■ . . j . / . / / ; 

20 • - • Suivant un exemple de realisation, le calcul numerique de l7mage d'uh 

cet operateur est effectue par exemple par resolution numerique de Tequatipn des oiides a 

. . ID pour le modele considere, en selectionnant des valeurs prises par le deplacement des 
paiticules aux emplacements, de capteurs et aux points d'echantillonnage en temps 
prealablement specifies et en appliquant un operateur propre a compenser les effets de 

25 divergence spheriqueet d' attenuation. - 

Suivant un -exemple de realisation, Toperateur numerique de modelisation du bruit 
est un schema numerique centre aux differences fin ies pour discretiser Tequation de 
transport du bniit et la fonction generatrice de bruit intervenant comme condition initiale 
en bordure de la dite zone associe a chaque bruit coixele, un espace. de fonctions 
30 generatrices du bmit foime par les fonctions du temps a support dans un intervalle de 
temps donne. 


' ' ' ' Differents exemples de normes ou semi-normes choisies pour Tespace des donnees 

et I'espace des fonctions generatrices de bruit seront donnes dans la description detaillee 
qui va suivre. 

Suivant un deuxieme mode de mise en oeuvre on recherche la distribution de 
5 perturbations, par rapport a un niilieu de reference prealablement choisi, de Timpedance et 
de la Vitesse dans la dite zone du milieu, les bmits correles affectant les donnees sont dus a 
des reflexions multiples dont la cmematique et les variations d' amplitude avec Toffset ont 
ete prealablement estiraes, les donnees mesurees sont captees par des capteurs sismiques 
de surface, I'emplacement des dits capteurs, le mode d'excitation sismique, la duree 
10 d'enregistrement et les points d'echantillonnage en temps etant definis, et Toperateur de 
modelisation etant defini via une linearisation de Tequation des ondes autour du milieu de 
reference. . 

Comme fonction cout quantifiant Tecart entre les mesures d'une part et la 
superposition de la reponse du modele et des bruits coneles associes aux fonctions 
15 . generatrices de bruit d' autre part, on chqisit par exemple le carre de la semi-norme de cet 
ecart dans Tespace des donnees. 

Suivant un exemple de realisation, I'ajustement du modele et des fonctions 
: . gen6ratrices du bmit est obtenu par une methode de relaxation par bloc pour eliminer les 
, inconnues correspondant a chaque fonction generatrice de bruits correles, cette methode de 
20 relaxation etant mise en oeuvre a Tinterieur des iterations d'un algorithme de gradient 
conjugue pour le calcul du modele. 

Presentation succiNCTE DES FIGURES 

Les caracteristiques et avantages de la methode selon Tinvention, apparaitront plus 
clairement a la lecture de la description ci-apres d'un exemple non limitatif de realisation, 
25 en sereferant aux dessins annexes ou : 

la figure 1 montre un exemple de donnees obtenues par mise eh oeuvre d'une methode 
de prospection sismique verticale (PSV) ; 

la figure 2 montre un modele de distribution varia tion d'impedance acoustique avec la 
profondeur ; 
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- la figure 3 montre des donnees de PSV .synthetiques obtenues sur la base du modele 
d'impedance de la Fig.2 ; 

la figure 4 montre un exemple de donnees de PSV contaminees par un seul bruit 
correle ; 

5 "la figure 5 montre un exemple de donnees de PSV contaminees par deux bruits 
coixeles ; ] . r , ; . " . 

^ - la figure 6 montre la distribution d'impedance en fonction de la profohdeur, obtenue 
, par inversion des donnees bruitees de la Fig.4 ; - , . 

- la figure 7 montre les residus d'inversion. (ecarts entre les donnees de^ la Fig. 4, et la.. 
10. reponse sismique du modele de la Fig ^6) ; , - . \ - . . ■ . 

- la figure 8 montre la distiibution d'impedance en fonction de la profondeur, obtenue^ 
par iriverision des donnees bruitees de la Fig.5 V ' , ' ; T 

- la figure 9 montre les residus d-inyersion;coiTespondants (ec^^ entre les donnees !de la; 
Fig. 5 et la reponse sismique du mbd^^ . '■ 1 

15 - la figure . 10 montre -la; distribution d'impedance en Tondti on de la profondeuir, obtenue ; 
■ ^ par inversion des donnees bruitees de la Fig;4 &n rech^rchanf le bruit correle sous la. 
foniae de la superposition de deux bruits coireles ayarit: des propnetes de propagation 
\ : ' ^ erronees c'est a dire differentes de celles du bruit apparaissant sur les donnees de la 

/pig- 4; • . ; ■ ■ • . . / : ^ 

20 r la figure 11 montre les residus d'inversion con*esp6ndants d'inversion (ecarts entre les 
^ donnees de la Fig. 4 et la reponse sismique du modele de la Fig. 10) ; 

- * les figures 12A, 12B montrent respectivement des donnees sismiques avec reflexions 

multiples, et la reponse sismique du modele obtenu apres inversion linearisee 
classique ; . : 

25 - les figures 13A, 13B montrent respectivement uri exemple de distribution d'impedance 

et de Vitesse de propagation auquel les donnees sismiques de la FigJ2A constituent la 
reponse sismique ; et 


' - les figures 14A, 14B montrent respectiveinent la comparaison des reponses sismiques 
des modeles obtenus par inversion classique (14A) el par inversion suivant la technique 
proposee incluant un pointe du "moveout" du multiple et une estimation des variations 
d'amplitude avec Toffset (14B). 

5 DESCRIPTION DETAILLEE DE LA ME THODE 

La formulation du probleme 

On dispose de donnees qui proviennent d'une mesure echantillonnee d'une 
fonciion. Cette fonction depend de plusieurs variables (variables d'espace ou d'espace- 
temps par example). 

10 Ces donnees, notees d, correspondent a des mesures effectuees pour nous 

renseigner sur un modele m. EUes se trouvent entachees de differents bruits : 
- des bruits correles (ou une superposition de bmits correles) ; 
des bruits non correles. 

Le probleme est de determiner quantitativenient le modele m (ou des fonctions de 
15 ce modele) a partir de la donnee d. 

Pour ce faire nous choisissons un operateur de modelisation F (lineaire ou non) qui, 
au modele m, associe la reponse du modele. Get operateur ne modelise en fait 
qu'imparfaitement le phenomene physique reel. C'est d'ailleurs essentiellement (mais pas 
exclusivement) pour cette raison qu'apparaissent, dans les donnees, des bruits con'eles : ces 
20 bruits correspondent a un signal lie au modele mais la relation entre Tun et I'autre nous 
apparait trop complexe pour etre incluse dans Toperateur de modelisation F que, pour des 
.. raisons diverses, nous voulons garder relativement simple. 

Dans ces conditions, nous modelisons le bruit en introduisant un . ou plusieurs 
espaces generateurs de bruit Bj (j variant de 1 a J) et les operateurs de modelisation de 
25 bruits associes que nous appelons Tj, Pour retrouver le modele m a pajtir des donnees d, 
nous proposons de rechercher ce modele comme solution du probleme d'optiinjsalion. 

min c{m,p,,p^,,.,i3j) 
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ou II 1^ est une norme (ou e vent uel lemon t une semi-nomie) dans Tespace des donnees. 

Le formalisme ci-dessus se distingue de Tapproche proposee par Nemeth et al 
( 1 999) de par le fait que: 

Toperateur de modelisation F n'est pas necessairement lineaire ; 

- nous Tious autorisons a considerer une superposition de bruits correles de types 
. differents (en ce sens qu'ils sont naodelises par des operateurs differents references par 

findice j), ceci au.prix d'une inflation des calculs tant en ce qui concerne le nombre 

- d'inconnues (ce qui complique singulierement la resdlution du probleme 
d' optimisation) qu'eri ce qui concenie le nombre de modelisations de bruit a effectuer ; 

- nous nous reservons la possibilite de chbisir a notre conyenance la norme | |^ ou 
eventuellemeW la semi -norme^ dans Tespace des donnees. ^ 

Le choix que nous allons faire pour | ||^ va .etre dicte pair des considerations liees^a 

refficacite de la methode de resolution (decrite plus bas), cette efficacite pennettant de' 
traiter le cas d'une superposition de bruits correles. II s'agit la d'une ouverture importante * 
6utre la possibilite de- traiter des dpiinees entachees d^ bruits correles alix caracteristiques 
complexes, on peut par ce biaiis surraonter les difficultes, rencontrees dans Tapproche de 
Nemeth et al,, liees a la^grande sensibilite du resultat. a une imprecision introduite lors de la 
definition; de roper^ de modelisation du bruit : il suffira, pour sunnonter cette 
difficult6, de rechercher le bruit par superposition de bruits associes a des operateurs/de 

modelisation Tj approximatifs. / 

' '.I " • . , . . . • I , - ■ • >' . f. 

' . . Notons que le choix du carre de la norme jj ||^ intervenant dans le probleme 
' d' optimisation n' a lien d'essentiel : nous ne changeons pas la solution en choisissant a la 
place une fonction qui, comme la fonction carre, est une fonction croissante sur 9\"^ (ou une 
fonction decroissante a condition de changermin en raa^^ ■/ 

La methode de resolution , 

La methode consiste a choisir convenablemcnt : 
la norme (ou semi -norme) || ||^ dans Tespace des donnees ; 

des normes || jj '. dans chaque espace Bj. ■ ' 


, de fa9on a. permettre la mise en place d'une solution algorithmique tres performante au 
probleme d' optimisation (1), Par "choisir convenablenient" nous en tendons faire en sorte 
que chaque operateur Tj constitue une isometrie (ou une approximation d'isonrietrie) pour 
les normes || ||^ et || ||^ respectivement dans I'espace de depart et dans I'espace d'arrivee. 

5 L'ldee directrice dans la determination de noimes permettant de rendre chaque operateur Tj 
isometrique repose sur I'existence d'egalites de type a conservation d'energie verifiees par 
les solutions des equations aux derivees partielles (ou d'egalites d'energie discretes pour 
les solutions de certains schemas nunieriques utilises pour la discretisation de ces 
equations). 

10 Si nous sommes capables de faire de tels choix, il est alors possible de rendre tres 

simple la minimisation de la fonction cout en eliminant les inconnues associees aux 
fonctions generatrices de bruits correles (determination du complement de Schur). Cette 
elimination est realisee a Taide d'une algorithmique adaptee telle que la methode de 
relaxation pai' blocs, un bloc etant associe a une fonction generatrice de bruits, ou la 

15 methode du gradient conjugue avec un preconditionnement de Gauss Seidel par bloc 
symetrique, etc., methodes de mise en oeuvre algorithmique toutes bien connues des gens 
de Tart, et que Ton trouvera presentees pai*: 

- . Golub G., H^et al (1983); "Resolution numerique des grands systemes lineaires" 

(EyroUes) 

20 La minimisation simplifiee de la fonction cout permet de reduire considerablement 

le temps necessaire a la resolution numerique. 

Premier exemple de mise en (Buvre applique a rinversion a ID de donnees PSV 
contaminees par une onde de tube 

Les donnees PSV (Pfbfil Sismique Vertical) font partie des acquisitions classiques 
25 de sismique de puits. Elles ont des applications multiples, Tune des principales etant leur 
utilisation pour faire le lien entre donnees sismiques de surface et mesures diagraphiques. 
L' inversion de donnees PSV telle que decrite par exemple par : 

- Mace, D, et Lailly, P. (1984) A Solution of an Inverse problem with ID Wave 
Equation applied to the Inversion of Vertical Seismic Profile; Proc. of the 6^^ Intern. 

30 conf. m "Analysis and Optimisation of Systems";'Nice, 2, 309-323 ; 
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a ete developpee a cette fin. Dans le probleme inverse ID, on suppose le modele de 
subsurface invariant lateralement et une excitation par ondes planes se propageant 
verticalement. Les inconnues du problenne sent la distribution d'impedance acoustique en 
fonction de la profpndeur (mesuree en temps de parcours vertical et sur un interyalle 
5 demarrant^ la profondeur du premier recepteur) ainsi que le mode d' excitation sismique 
(de fagon precise la fonction du temps caracterisant la condition aux li mites a la profondeur 
du prernier capteur). II s'agit la d'un problerne inverse non lineaire, les donnees PSV 
dependant de fagon pon lineaire de la distribution d'impedance. 

Les donnees PSV sont souyent contaminees par un bniit cdrrele de tres forte 
10 arnplitude : Tonde de boue. II s'agit d'une onde guidee dans la boue envahissant le puits. 
Sa Vitesse de propagation est lente par rappprt a la vitesse de propagation dans les roches a 
. I'entour dii puits. De plus la vitesse de propagation de cette onde depend de la frequence 
(propagation dispersive). Enfiii cette onde se reflechit notamment\au fond du puits, 
donnant de plus naissance a un autre mode de propagation. Des donnees PSV typiques spnt 
15 . donnees sur la figure ci-dessous :' on y voit les principal e? caracteristiques de I'onde de'; 
tube/ Son amplitudp et Textension de la zone contaipinee rerident difficilement exploitable^ 
ie signal contenu dans les nombreux echantillons sismiques cqntamines par Tonde de.boue. I 

Nous decrivoris ci-apres la fnise en ceuvi-e; de k methbde.-pCur r inversion a ID des . 
'■' "donnees PSV. - ^ ■ , ' ^ y ■ 

20 Le p 'rdtqcole experimental 

II consiste en la specification : 
- ' des differentes cotes auxquelles les recepteiirs seront places dans le puits (suppose 
, vertical) : ;dans nos experiences ies capteurs couvrent i'intervalle de profondeurs 
mesurees en temps simple [Xmin= 0.05s ; Xmax= 0:2s] , avec un recepteur tous 
25 les Ax=lms (en temps simple) ; nous.obtenons airisi de echantillons, numerotes de 0 a 

de la nature physique de la mesure .effqctuee par ces capteurs : dans nos experiences les 
capteurs mesurent le deplacement vertical resultant de la propagation des ondes ; 
de la duree des enregistren^ents : dans nos experiences les capteurs mesurent Tetat 
30 vibratoire sur I'intervalle de temps [0 ,T=1^5s], ces donnees etant echantillonnees tous 

les At=lms. Nous obtenons ainsi des echantillons numerotes de 0 a N. 


Nous appelons dj" rechantillon de donnees enregistre a la profondeur Xniin+ i Ax et 
au temps n At. Nous avons bien sur : T= N At et Xniax= Xjnin+I Ax. 

Acquisition des donnees conformes a ce prowcole experimental 

Nos experiences ont ete faites a partir de donnees synthetiques c'est-a-dire obtenues 
5 par modelisation numerique, elle meme reaJisee en deux etapes : 

Creation de donnees synthetiques non bruitees : celle-ci a ete realisee par resolution 
numerique de Tequation des ondes acoustiques a ID, le modele etant la distribution 
d'impedance acoustique representee sur la Fig. 2 (la profondeur est mesuree en temps 
simple) et le mode d' excitation (condition aux limites de Neumann a la profondeur 0) etant 
10 une classique ondelette de Ricker. 

. Les donnees PS V synthetiques issues de cette modelisation sont representees dans 
. la figure 3 : I'axe vertical representant le temps d' observation est gradue de 0 a 1,5 s et 
r axe, horizontal representant la profondeur (en temps de parcours simple) des differents 
capteurs est gradue de 0,05 a 0.2 s. 

15 Adjonction de bruit sur les donnees 

Nous allons perturber les donnees PSV ainsi calculees en leur ajoutant d'une part 
un bruit aleatoire (mais filtre dans la bande des frequences sismiques) d' amplitude assez 
forte et d'autre part un ou plusieurs biuits coireles d' amplitude tres forte. Chacun de ces 
bmits se propagera vers le bas avec une vitesse lente qui depend de la frequence 

20 (propagation dispersive) : ces bruits correles sont censes representer des ondes de boue. lis 
ont ete modelises par resolution numerique (utilisation du schema centre aux differences 
finies) d'une equation de transport avec une vitesse de propagation constants La 
dispersion des ondes a ete obtenue en utilisant dans de grands pas de discretisation dans le 
schema aux differences. Dans le cas de la superposition de plusieurs bruits coiieles, les 

25 ; vitesses de propagation associees a chacun de ces bruits sont differentes. Les deux figures 
4, 5 montrent les donnees PSV ,contaminees par un seul bruit conele. (Fig.4) et par la 
superposition de deux bruits correles (Fig.5). 
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Specification de Voperateur de modelisation 

La modelisation est realisee par la resolution de T equation des ondes ; 
a\TBC ia coTidition 'anx limites de Neiimann : 

.. . . • . . , . . * ■ • • ()y ' . ■■ - 

. ■ ■■ . ' ' • : -rr(^min)^(j:mm, *) = h{t) sur [0, T] 

. et les conditions initiales : - 

' La definition de rbperaiteur de modelisation demande la definition prealable cle 
5 ropefateur d'observation qui fom£^^^^ L'operateur' 
d'observation est done roperateur qui; .a lai fonGtion y(x,t) solutioii du systeme ci-dessus, 
associe les echantillons y(xut") p:o 

* ■ L'operateur de mbdelisatiofi que nous appellerons F(m) est done roperateur'(non' 
. Jineaire) qui au couple de fonctioris' in=(a(x)Ji(t)) assod echantillons yCxj^t") pour; 
10 /; Xi=Xmin+i Ax et t"-h At, i^^O,.:^ ■ / ^ ' 

, ' Specification de ropemteur de modelisation associe 

\ - Nous precisons dans ce paragraphe la procedure pouvarit etre \utilisee pour 

' ' ■ modeliser chaque bruit correle. Chaqiie briiit correle sera caracterise par des directions de 
correlation qui lui seront specifiques et que nous supposoiis connues; au moins 
15 . : approximativement : ces directions de correlation sont.specifiees par rintermediaire d'un 

; champ de vecteurs de coirelation cj (x,0^de conriposantes (cj'^CXjt), CjVx,t)) qui devra etre 

defini en tout point du domaine I^„^„,Jc„u,x]^[f^'^]- Nous pouvons de fa9on equiyalente 

specifier les directions de con;elation par Tintennediaire d'un faisceau de lignes de 
correlation qui representera les lignes de champ associees a c.(xj). De telles lignes de 

20 correlation peuvent etre obtenues suivant la technique decrite dans les brevets EP 354.112 
(US 4 972 383) du demandeur, a partir d'un pointe de quelques phases caractefistiques du 
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bruit rinformation etant etendue a tout le domaine Un^io^-^maxl^fo^r] en utihsant des 

procedures classiques d 'interpolation et/ou d'extrapolation. La specification des directions 
de coiTelation revient a specifier la distribution de vitesse de propagation du bruit Cj(x,t) 

c'^fx t) / 

qui est liee au vecteur de coirelation par la relation c.(x,t)= ^ Nous 

supposons ici qu'il n'y a pas de variations d'amplitude de Tonde lors de sa propagation 
(une situation differente sera presentee pour le 2^'^^ exemple d' application). La 
modelisation d'un bruit correle tel que specific ci-dessus est basee sur la remaique 
suivante : une onde se propageant le long des lignes de coirelation sans changer 
d'amplitude est solution de I'equation de transport : 


^b{x, t) , c{x, t) = 0 


Dans ce contexte, nous introduisons 

- un espace Bj de fonctions generatrices de bruit qui sera Tespace des fonctions j3j(t) a 
support dans r ]cz [o,r] et que nous prolongeons par 0 sur I'intervalle [0,r] 
entier ; 

.15 - Toperateur de modelisation du bruit Tj : 

Tj-: fijit) e Bj ^ bj{x, t)eD 
solution de I'equation de transport 

■ " ' ' et Y^rifiant les t!onditions initiales : 


^ hj{x = a:^.,n, f) = Bjif) 
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II faut remarquer que la geometrie du faisceau de lignes de correlation (ou, ce qui 
revient au meme, la distribution.de vitesse de propagation) ne peut pas etre quelconque : 
les lignes de correlation ne sont pas autorisees a se couper : sinon le probleme de resolution 
de Tequation de transport serait un probleme mal pose. Pour les niemes raisons les lignes 

5 de correlation ne sorit pas autorisees a intercepter deux fois les bords sur lesquels les 
conditions initiales sont specifiees. Gette demiere consideration peut amener (pour, 
inodeliser un bruit correle coiTespondant a une pnde montante), a specifier la condition 
initiale non plus en x-Xnun mais -en x=Xmax on merne en une valeur interriiediaire quitte a 
decomposer le probleme de Cauchy en deux problemes associes aux domaines situes de 

iO part et d'autre de ce bord inteiroediaire, ' . . 

Nous supposons par ailleurs que la geometrie des lignes de correlation est telle' que 
(des situations plus generales peu vent etre vCnyisagees mais les resultats qui suivent 
. prennent aiors une fontie dif^^^ . 

. ' ' ■ ■• ' ' ^ 

15 " ;Enfin ''nous supposons que le faisceau de ligneV de:0r^^^ qui intefceptent en ; 
. ; x=Xmin'l'intervalle;{r j""'" , r^"^' ] n'interceptent pas en t~T r-intervalle [x^^^ , x^^ ] . . 

En fait Fequation de transport est resolue par rinternaediaire d'un schema 
numerique. On peut par exemple utiliser le. classique schema centre aux differences finies. 
^Pour ce faire nous choisissons des. pas de discretisation Axj' et Atj* (qui ne sont pas 

20 necessairepient egaux a Ax et At mais que nous supposerons, meme si des situations plus 
generales peuvent $tre envisagees, sous multiples de ces quantites) et introduisons un 
maillage dont les ncxuds sont les points de xodrdpnnees (Xniin+ i''Axj',n'Atj') avec i' et 
neN, Le classique. schema centre aux differences finies explicite ci-dessous permet, 
partant des conditions initiales, le calcul de prpche eh prpche des differentes valeurs prises 

25 par la fonction bj(x,t) aux differents noeuds du maillage pour sirnplifier les notations nous 
avons omis I'indice j associe au bruit coirele considere). 
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Dans la formule ci-dessus la quantite Oi-^ represente revaluation de la quantite a 
(a est un terme generique) au point de coordonnees (XniiTi+ i'Ax,j',n'Atj') 

Bien entendu Tutilisation de pas de discretisation Axj' et Atj' suffisamment petits 
5 devant la longueur d'onde est requise si nous voulons que le schema numerique donne une 
approximation precise de la fonction solution de Tequation de transport. 

Cependant, I'utilisation de pas de discretisation plus grands (mais toujours plus 
petits que la longueur d'onde) permet de modeliser des phenomenes propagatifs plus 
complexes en rendant la vitesse de propagation de Tonde dependante de la frequence : on 
10 peut ainsi modeliser des propagations dispersives telles la propagation des ondes de tube. 
Le choix de pas de discretisation adaptes a la modelisation d'un mode particulier de I'onde 
de tube peut etre realise : 

En determinant, par un examen des donnees, la vitesse de propagation de ce mode 
(qui depend ici non seulement de I'espace et du temps mais encore de la frequence) ; des 
15 techniques de tomographic peuvent etre utilisees a cette fin, telles que celles presentees 
dans la publication suivante: 

- Ernst F. et al (2000); Tomographt of Dispersive Media; J. Acoustic. Soc. Am., 108, 
105-116. 

En connaissant les proprietes de propagation dispersive associees au schema 
20 numerique, proprietes qui decoulent d'une analyse de la relation de dispersion numerique 
(comme decrit pai' exemple dans les publications suivantes: 

- . Alford R. M. et al, (1974), Accurency of Finite Difference Modeling of the Acoustic 

Wave Equation; Geophysics, 6, 834-842. 


1=^ 
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Specification d'une norme ou d'une semi-norme dans Vespace des donnees 

Un premier choix consiste a prendre comme nonne dans Tespace des donnees une 
norme discrete quelconque qui se veut une approximation (via une formule de quadrature) 
de la nomie defiiaie en continu par : 


1 


• Nous pouyons aussi, ce qui cdnstitue, cornme notis le veaons.au § suivapt, un 
meilleur choix, utiliser la seriu-nprme : : ' ' - ' . . ' ' ' 


II 


. . -A- . T==0 .77=0 ^ _ r 


(2) 


10 . oij u est un yecteur quelconque de I'espace <des^^^^ 

. . les nunieros d'echantiUons respectiveni ' . ^ : / 

Pai- ailleurs d'autres choix sent possibles. ^ ' ; 

^ . ■ Specification, dans chaque espace Bj de fonctions generatrices du bruit, d'line 
■ nonne pour laquelle chaque operateur de niodelisation du bruit Tj 'etablit/une relation 
15 isometrique (ou V approximation d'une relation isometrique) entre Vespace des fonctions 
generatrices du bruit et Vespace des donnees (muni de la semi-nonne definie en 5) 

Nous precisbns dans ce paragraphe la procedure pouvant etre utilisee pour definir la 
noiTiie dans ] 'espace des fonctions generatrices associe a chaque bruit coiTcle. Cette 
procedure etant la meme pour chaque bruit, nous la decrivons de fagon generique en 
20 omettant, dans les formules, Tindice j associe au bruit considere. 
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Un premier choix consiste a munir chaque espace B de foncdons generatrices de 
bruit d'une noiTOc discrete quelconque qui se voudrait une approximation (via une formule 
de quadrature) de la nonne definie en continu par : 


(t) 


5 auquel cas, en utilisant le fait que, grace aux hypotheses enoncees au §4, on a bj(x,T)=0 
pour ^' ^ [^niin ' -^max ] ' Toperateur lineaire Tj realise une isometrie entre Bj et D. Mais, 
comme le calcul de Tj/3j est effectue numeriquement. le caractere isometrique de 

I'operateur T ne sera pas parfait mais seulement approche, Tapproximation etant d'autant 
meilleure que les pas de discretisation Aa% A/, Ax^ \Atj 'sdnt petits. Le fait que nous n'ayons 

10 alors qu'approximativement une isometrie se traduira par des performances moins bonnes 
au niveau de la mise en ceuvre algorithnfuque presentee au paragraphe 8 ci-apres. 

Un meilleur choix consiste a munir chaque espace Bj de fonctions generatrices de 
bruit de la semi nonne : 


m\B - A.:A^ / ~r {p-"^' + rf 

(3) 

15 Dans ce cas I'operateur Tj realise rigoureusement une isometrie entre Bj et D, du 

moins si = 0 Vz = 0,/ - 1 . Cela resulte de Tegalite d'energie discrete : 

. n=0 ' i=Q n=0 . " 

^gaJit^ valable dans le cas (sirnplifi^ pour les besoius.de la presentation) on 
Ax' Ax et At' = At et si Fon a d^fini : f i — f i 
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Definition de lafonction cout 

Nous definissons la fonction cout par la formule : 

it 2 

Recherche du mddele et des fonctions generatrices du bruit minimisant lafonction 
' 5 '■■ cout. cette recherche se faisant par une methode qlgdrithmique tirant parti, explicitement 
ou implicitement, des proprietes d' isbmetrie des operateurs de modelisation dU bruit (cf6). 

. 'La methode algoiithmique exploite la speciflcite du complement de Schur, 

specificite assbdiee au caractere isometriqu^ des operateurs Tj, Ceci peut .se faire par 
; exemple en. realisant que . mininiisdr revient & . . minimjser 

10 C\m) = C(m;^i(?^),y52(m),...J; (m)) on \ ' (^^^ ■ minimis^e^ 

C{m,pyyfi^^...pj) ponx m dorine. Notdns :que la hessienue associep a cette forme 
quadratique est le complement de JSchur. La minimisation de C\m) peut etre realisee par 
■ : une quelconque ';meth6db; d'optimisatipfi telle ^^^^^^ BFGS d' une. methode dei 

/ . Quasi Ne>vton, .-, V'^ V' ; ' ' : • - >: 

15 La deterniination de (y§j(m),y52(^)'-^^^^ aisee grace au caractere 

isometrique des operateurs Tj., notaminent-si nous avons effectues les choix (2) et (3) 
respectivement pour definir |||^ et . S'il.n'y a qu'un seul type de bruit (J=l) et si 

\ I'isometrie est parfaite, la determination de ^^(m) demapde simplement une evaluation du 
■gradient de la forme quadratique. S'il y a superposition de plusieurs bruits (J>1), differents 

20 algorithmes peuvent eti'e envisages pour profiter du caractere isometrique des operateurs 

Tj. On peut par exemple detenfiiner {Pi^{m),p^{m),...^ j{rn)) par une methode de 

relaxation (alias Gauss-Seidel) par bloc, un bloc etant associe a un groupe crinconnues 

pj{m) , une iteration de relaxation ne demandarit simplement (toujours dans le cas d'une 

iisoiTietrie parfaite) qu'une seule evaluation du gradient de la forme quadratique associee au 
25 bloc considere. Pour les problemes rendus difficiles par la proximite des directions de 

con'elation des differents types de biTiits/ on peut envisager une methode de gradient 
conjugue preconditionne, la matrice de preconditionnernent etant de type Gauss-Seidel par 


zu 

bloc (alias relaxation) symetrique. La encore la mise en oeuvre du preconditionnement est 
tres facile, la resolution du probleme associe a un bloc ne demandant qu'une evaluation du 
gradient de la forme quadratique associee au bloc considere. 

Les Fig.6 all montrent les resultats obtenus sur trois types d'experiences de mise 
5 en oeuvre de la methode. 

Dans le cas dWe inversion de donnees bixiitees (Fig.4) a la fois par la super-position 
d'un bruit correle et d'un bruit aleatoire, la Fig.6 donne la distribution d'impedance 
retrouvee, tandis que la Fig.7 montre les residus d'inversion. 

Dans le cas d'une inversion de donnees bruitees (Fig.5) a la fois par la superposition 
10 de deux bruits correles et d'un biiiit aleatoire, la Fig.8 donne la distribution d'impedance 
retrouvee, tandis que la Fig.9 montre les residus d'inversion. 

Dans le cas d'une inversion de donnees bruitees (Fig.4) a la fois par la superposition 
d'un bruit correle et d'un bruit aleatoire ou Ton ne connait pas de fagon precise les 
proprietes de propagation de I'onde de boue, il est nature! de rechercher un modele de bruit 
15 correle constitue de la superposition de deux bruits coixeles ayant des viitesse de 
propagation erronees mais encadrant la vraie vitesse de propagation de I'onde de boue. La 
Fig. 10 donne la distribution d'impedance retrouvee, alors que la Fig, 11 montre les residus 
d'inversion. 

Deuxieme exemple d*appIication de la methode a la determination d'un modele par 
20 inversion linearisee de donnees multi-deports (multi-offsets) contaminees par des 
reflexions multiples 

Cette application se distingue de la premiere principalement par le fait qu'elle 
illustre les possibilites qu'offre la methode vis a vis de la prise en coinpte de variations 
d'amplitude du bruit le long des directions de correlation. Une.autre difference reside dans 
25 les proprietes de propagation des bruits correJes, la propagation des. reflexions multiples 
n'etant pas dispersive. 

L'inversion linearisee figure parmi les methodes sophistiquees utilisees par les 
geophysiciens pour obtenir un modele quantitatif des heterogeneites de la subsurface, ce 
type d' information etant tres precieux notamment pour la caracterisation des reservoirs. 
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Les donnees sont ici des donnees sismiques de surface. Cette methode est decrite par 
exemple dans la publication suivante : . 

Bourgeois A. et aL, (1989) "The Linearized Seismic Inverse Problem: an Attractive 
Method for a Sharp Descnptiopn of Strategic Traps": 59^^" Ann. Intemat. Mtg. Soc. 
5 Expl Geoph. Expanded Abstract, pages 973-976. 

■Par ailleurs il est necessaire- de disposer d'un modele de reference, constitue (dans 
■ le oontexte d*;un modele acousti que) d'une distribution de vitesse et d*une distribution 
d'impedance acoustique : ce sera le module autour duquel on effectuera une jinearisation 
dans I'equation des ondes (approximation de Bom). Pour simplifier la presentation nous 
10 ' nous placerpns dans le cas ID mais cette hypothese e$t d'autant moins essentielle que les . 
applications pratiques concement pour I'essentiel la sismique 3D. Dans le contexte'- ip, 
toute rinformation sismique. est contenue dans un seul point de tir: Par ailleurs le modele. 
de reference est lui aussi ID :t il ne depend que de la prpfondeiir. Les inconnues dii; 
- probleme sont la distribution de ^pertiirbation d-impeda;itQe acoustique et de vitesse en 
- . 15. fonctiqn jde la profpndeur. En revanche, contrairement au'pfemier exenriple d'appHcation le; 
mode d'excitation sismique (et notammeht 1 -dndelette sisraique) seront.sUpposes connus. Hi 
: s'agit la d'un pfrobleriie inverse rendu lineaire grace a la lineaiisation. , :, 

Les- donnees sismiques de surface ;S&nt: sou vent par des/reflexions/ 

inultiples qui peiivent etre d'assez forte amplitude. Leur cinematique, qui, dans le contexte 
20 ID, se baracterise par le "moveout", est en general assez differente de celle des reflexions 
primaires : ces deux types d'ondes ont, en general, traverse des couches geologiques ayant 
des vitesses de propagation differentes. Voila pourquoi Tin version Jinearisee constitue en 
. elle-meme, un technique efficace pour attenuer les reflexions rnultiples : la cinematique 
des reflexions primaires (les evenements modelises paivl' approximation de Bom) est 
25 entierenient definie par la distribution de vitesse dans le modele de reference: les 
evenements modelises ne pfeuvent done s'ajuster suf des reflexions multiples pour peu que 
le "mdveout" de ces demieres differe de celui des' primaires. Cependant, en raison des 
stationnarites a Toffset nul, il y a, en pratique un ajusterhent partiel et les resultats 
d'inversion lineaiisee s'averent contamines par la presence du bruit coiTcle que constituent 
30 les reflexions multiples. Les donnees sismiques de la Fig,12A ont ete modelisees avec 
prise en compte des reflexions multiples (il n'y a que 3 . primaires qui arriyent a Toffset nul 
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successivement aux temps :500 ; 1000 ; 1650 ms, la 2^"^^ arrivee etant d'amplitude 
particulierement forte). 

Sur la reponse sismique du modele obtenu apres inversion linearisee classique (cf. 
Fig.l2B), le multiple de forte amplitude arrivant vers 1500 ms a I'offset nul n'a ete que 
5 faiblement attenue et le modele se trouve Gontamine par le multiple, II est alors naturel de 
tenter d' exploiter les possibilites offertes par la methode de Nemeth pour mieux 
discriminer multiples et primaires, notamment avec les ameliorations apportees par notre 
methode. 

Nous decrivons ci-apres la mise en oeuvre de la methode pour Tin version linearisee 
10 a ID de donnees sismiques de surface. 

Le protocole experimental 

II consiste en la specification: 

Du mode d' excitation sismique : dans nos experiences la source sismique est un point 
source module en temps par une fonction (ondelette sismique) notee 'w(t) ; 

15 - des differents emplacements auxquels les recepteurs sont disposes : dans nos 
experiences les capteurs sont sitiies a la prof ondeur 10m et couvrent Tintervalle 
d'offsets [Xrtiin-0 ,Xmax=J500m], avec un recepteur tous les Ax= 100m; nous obtenons 
' ■ ainsi des traces sismiques numerotees de 0 a I ; 

- les capteurs mesurent le champ de pression resultant de la propagation des ondes ; 

.20 • - de la duree des enregistrements : dans nos experiences les capteurs 'mesurent I'etat 
vibratoire sur I'intervalle de temps [0 ,T=1 800ms], ces donnees etant echantillonnees 
tous les At=lms. Nous obtenons ainsi des echantillons numerotes de 0 a N. 

Nous appelons di" I'echantillon de donnees enregistre pour I'offset x,nm+ i Ax et au 
■ . temps n At. Nous avons bien sur : T= N At et Xjnax= Xmin+I Ax. 

25 Acquisition des donnees conformes a ce protocole expenmental 

Nos experiences out ete faites a partir de donnees synthetiques : ces donnees ont ete 
obtenues par resolution numerique (methode des differences finies) de ] 'equation des 
. ondes 2D : 
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Hxec la condition aiix limites : 

P{z^ ())==() 

et les conditions initiates : 

equation dans.laquelle : . • . . 

- X, z, 1 4esignent respecti vement la coordoniiee laterale, la profondeur et le temps ; 

- at c representent respecti vement les ^distributions Xfonctions de la profondeur) 
-5 . d'irnpedance acoustique et de Vitesse de propagation 

L'ondelette choisie est -une classique brideiette de Ricker de frequence centrale 
30Hz. et les distributions d'impedance a et de Vitesse c soht donnees sur les figure 13A, ■ 
13B. La reponse sismique de ce modele est done constituee des donnees representees sur la 
; Fig. 12A. Nous appelons Fnl I'operateur qui, au^couple de fonctions (a(z),c(z)), associe les : 
10 donnees sismiques. ' . . . ' c. 

Specification de V operateur de modelisati^^ , ^ - . ; 

■ ^ Nous choisissohs un modele de reference decrit par le couple de fonctions 
. (ao(z),cp(2)), pans les experiences presentees dans la suite nous avons choisi co(z)^c(z) et 
ao(z)=cte (=1). Nous choisissons pour operateur de modelisatipn I'operateur jacobien de 
15- I'operateur Fnl au point -(ao(z),cp(z)). En fait, pui^que nous avons choisi co(z)=:c(z), la 
seule inconnue est 5a(z)= a(z)- GoCz) et seule la composante corresporidante du jacobien 
est irhportante. ^ 

L'operateur de modelisation que nous appellerons F(5iTi) est done Toperateur 
(lineaire) qui au couple de fonctions 6m=(Sa(z),5c(z)) associe les echantillons 5y(Xi,t") qui 
20 sont les composantes du vecteur F (5m) pour Xi^Xmin+i Ax et t"=n At, i=0, . . . ,1 ; n~0, . . . ,N. 

Specification de Voperateiir de modelisation associe a chaque bruit correle 

Nous precisons dans ce paragraphe la procedure pouvant etre utilisee pour 
modeliser chaque bruit correle. Chaque bruit correle sera caracterise par des directions de 
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correlation qui lui seront specifiques et que nous supposons connues, au moins 
approximativement : ces directions de correlation sont specifiees par Tintermediaire d'un 
chanap de vecteurs de correlation Cj(xJ) de composantes (cj^Cx,!), Cj\x,t)) qui devra etre 

defini en tout point du domaine Uj^n'^maxl^lo^^]- Nous pouvons de fagon equivalente 
5 specifier les directions de correlation par Tintermediaire d'un faisceau de lignes de 
correlation qui representera les lignes de champ associees a Cj{x,t), De telles lignes de 

correlation peuvent etre obtenues a paitir d'un pointe de quelques phases caracteristiques 
du bruit rinformation etant etendue a tout le domaine U'^in, Jf^maxl-^lo*^] en utihsant des 
procedures classiques d'interpolation et/ou d'extrapolation, comme decrit dans le brevet 

10 . EP 354.112 du demandeur deja cite. Pour ce deuxieme exemple d* application dans lequel 
nous voulons eliminer le multiple precedemment mentionne, nous pointons ce multiple ( 
par exemple le pic d' amplitude), ce qui nous permet de definir les variations de temps 
d'arrivee en fonction de Toffset et nous definissons un continuum de lignes de correlation 
par une simple translation verticale de la ligne pointee (dans ces conditions le vecteur 

15 c(xj) ne depend pas de t). La specification des directions de correlation revient a specifier 
la distribution de vitesse de propagation du bruit Cj(x,t) qui est liee au vecteur de 

c^ (x t) / 

corr61ation par la relation c . (x,t) = ' / . , ^ 

' / c;(x,t) 

Contrairement au 1^^ exemple d'application deja deciit, nous avons ici de variations 
d'iumplitude du bruit le long de la ligne de con-elation. Nous supposons que nous pouvons 
20 estimer ces variations d* amplitude a partir d'une mesure sur les donnees ou bien a Taide de 
considerations theoriques : cette mesure nous defmi une fonction g(x). 

La modelisation du bruit se fera done par la composition de deux operateurs : 

Un operateur de transport, comme dans le premier exemple d'application ; 

Une modulation d'amplitude, c'est a dire la multiplication de la solution de 
25 I'equation de transport par la fonction g(x). 

En reprenant le schema donne pour le 1^^ exemple d' application (et notamment en 
faisant appel aux memes hypotheses et notations), nous introduisons : 
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un espace Bj de fonctions generatrices de bruit ciui sera Tespace des fonctions fijit) a 
support dans [r"'"\/;)"''']c [o,r] et que nous prolongeons par 0 sur I'intervalle [o,r] 
entier ; 

I'operateur de modelisation du bruit Tj : 

, V " solution de r^quation de transport. 
■ ■■ ■ • ^hj{x,t)%Cj{x.,t) = 0 dank[xir,.in. r^mar] X [0,T] 

et v6rifiant les conditions initiales : 


. t - 0) = 0 ^ . 


; En fait i'equation de transport est resoiue par il'in^^ d'un schema 

numerique. On peut par exemple litiliser le classique schema centre aux differences finies. 
Pour pe faire iious cho^ et Atj' (qui ne' sont pas; 

necessairement .egaux a Ax et At mais que iious supposerons, meme si des situations plus 
generales peuvent etre enyisagees, sous multiples de ces quantites) et introduisons /uii 
maillage dont les noeuds sent les points de coordonnees (XiTTjn+ i'Ax/,n'Atj') avec i'- et 
n'E A/' . Le classique schema centre aux differences finies explicite ci-dessous pemiet, 
paitant des conditions initiales, le calcul de proche en proche des differentes valears prises 
par la fonction bj(x,t) aux differents nceuds du maillage (pour simplifier les notations nous 
avoiis omis Tindice j associe au bruit correle considere). 


Dans la formule ci-dessus la quantite Qi^"' represente revaluation de la quantite a 
(a est un terme generique) au point de coordonnees (Xmin+ i' Axj',n' Atj'). 

Ici les pas de discretisation Axj' et Atj' doivent etre choisis suffisannnent petits 
devant la longueur d'onde car nous voulons que le schema numerique donne une 
approximation precise de la fonction solution de 1' equation de transport (nous ne voulons 
pas de dispersion). 

La demiere etape dans la modelisation du bruit consiste a selectionner les 
echantillons bj*" qui appaitiennent aux noeuds communs aux deux maillages et a les 
multiplier par g(iAx) : nous obtenons ainsi les differentes composantes du vecteur (de 
Fespace des donnees) Tj(j3j) , 

Specification d'une nonne ou d'une semi-norme dans Vespace des donnees 

Un premier choix consisterait a prendre comme norme dans I'espace des donnees 
une norme discrete quelconque qui se voudrait une approximation (via une formule de 
quadrature) de la norme definie en continu par : 



Nous pouvons aussi, ce qui constitue, comme nous le veiTons au § suivant, un 
meilleur choix, utiliser la semi-nome : 

1 

■ ■ ■ \ ^ 1=0 n=0 ^ ' ■ ' I ' 

/ ' ' " (4) 

oii u est un vecteur quelconque de I'espace des donnees (i et n sont les indices representant 
les numeros d' echantillons respectivement en espace et en temps). 

Par ailleurs d'autres choix sont possibles. 
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Specification, dans chaque espace de fonctions generatrices du bruit, d'une norme 
pour laquelle chaque operateur de modelisation du bruit etablit une relation isometrique 
(on V approximation d'une relation isometrique) entre V espace des fonctions generatrices 
du bruit et V espace des donnees (muni de la semi-norme definie au § 5) 

5 Nous precisons dans ce paragraphe la procedure pouvant etre utilisee pour definir la 

norme dans 1' espace des fonctions generatrices asspcie a chaque bruit correle., . 

Un premier choix consisterait a munir chaque espace Bj de- fonctions generjatripes 
de bruit d'une norme discrete quelconqjue qui se voudrait une approximation (via une 
. formule' de quadrature) jit la norme definie eh cpntinu par (la encore ripus aVQns omis 
10 I'indice j pour aI16ger les notatipns) : . ... V, 


\\in\B = ^ 


auquel cas, en utilisant le fait .que bj(x,T)=0 pour;>'€ [^niiu , ^^ax ] ' 1' operateur lineaire Tj \ 

realise une isometrie enfre Bj et /D. Mais, cdmrne; le calcul dc T^jj3j est . effectu^ . 

numeriquement ' le cai^actere isometrique de l^operiateur Tj ne . sera parfait' mais 
15 seulement approche, r approximation etant d'autant meillcure que les /pas de 
discretisation Ax, A?, Ax^ A? ^ 'sont petits. Le , fait que nous n^ayons alors 

qu'approximativement une isometrie se traduira par des performances moins boiines au 
niveau de la mise en oeuvre algdrithmique presentee au paragraphe 8. - - 

Un meilleur choix consiste a munjr chaque espace Bj de fonctions generatrices de 
20 bruit de la senoi norme (la encore nous avons omis rindice j pour alleger les notations) : 


/-I \ /;V-1 


\W\\b = . AxSd [Y^gHiAx)] lY^^i^n+i^^n 

\ \i=0 / \n=Q ^ ~ 


(5) 
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Dans ce cas I'operateur Tj realise rigoureusemenl une isometrie entre Bj et D, du 
moins si la solution dii schema de discretisation de requation de transport est 
identiquement nulle pour n=N, 

Definition de lafonction cout 

5 Nous definissons la fonction cout par la foimule : 

Notons que, dans la formule ci-dessus, FhlOtcio) est un vecteur constant. 

Recherche du modele et des fonctions generatrices du bruit minirnisant lafonction 
cout, cette recherche se faisant par une methode algorithmique tirant parti, explicitement 
10 ou implicitement, des proprietes d' isometrie des operateurs de modelisation du bruit 

(cfm 

La methode algorithmique exploite la specificite du complement de Schur , 
specificite associee au caractere isometrique des operateurs Tj. Ceci peut se. faire par 
exemple en realisant que minimiser C{Sn, , y?^ ,„.pj ) revient a minimiser 

15 C\Sm) - CiSmJ,(dm)J^{SmX...^j {&n)) ou 

ifi^iSm), P^{5m),,,.^j (Sni)) minimise C{dm, pour 5m donne. A noter que la 

hessienne associee a cette forme quadratique est le complement de Schur. La minimisation 
de C\Sni) peut etre realisee par une quelconque methode d'optimisation ; la fonction cout 
C\5m) etant quadratique, I'utilisation d'une methode de gradient conjugue est ici tout a 
20 fait indiquee. 

La determination de i^i(m)//3^(m),,„^j{m)) est tres aisee grace au caractere 
isometrique des operateurs Tj,, notamment si nous avons effectues les choix (4) et (5) 
respect] vement pour definir 1 1|^ et || ||^ , S'il n'y a qu'un seui type de bruit (J=l), ce qui est 

le cas pour notre illustration el si risometrie est parfaite (ce qui est egalement le cas grace 
25 aux choix que nous avons effectues) la deteiinination de fi^{m) demande simplement une 
evaluation du gradient de la fonne quadratique. S'il y a superposition de plusieurs bruits 
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(J>1), differents algorithmes peuvent etre envisages pour profiler du caractere isometrique 
des operateurs Tj. comme expose pour le premier ex emple d' application. 

La Fig. 14B montreftt les resultats obtenus (reponse sismique du modele solution) 
par mise en oeuvre de la methode-dans le cas ou la fonction g(x) est choisie sous forme tres 
5 simple (fonction affine). On peut constater T amelioration par rapport au resultat 
d'inversion classique (Fig. 14A).. 

On a decrit .des exenntples de niise en oeuvre bu le parametre physique dont on 
modelise la distribution dans le spuS'Sol est I'impedance acoustique. D est bien evident que 
Ja riietHode dans sa definition la plus generale, peut s'appliquer a la recherche de la 
10 distribution de grandeurs physiques affectes par des bruits correles dans tout milie.u 
heterogene. 


REVENDICATIONS 

1) Methode pour estiiner, a partir de donnees obtenues par exploration d'une zone 
d'un milieu heterogene, un modele representatif de la distribution, dans la zone, d'au moms 
une grandeur physique, affranchi de la presence de bruits con-eles pouvant entacher les 

5 donnees, caracterisee en ce qu'elle comporte les etapes suivantes : 

a) acquisition de mesures renseignant sur certaines caracteristiques physiques de la 
zone en suivant un protocole experimental defini ; 

b) specification d'un operateur de modelisation, qui, a un modele de chaque 
grandeur physique associe des donnees synthetiques constituant la reponse du modele, les 

10 mesures et les donnees synthetiques appaitenant a I'espace de donnees ; 

c) pour chaque bruit coirele, repere par un indice j variant de 1 a J, choix d'un 
operateur de modelisation qui associe un bruit correle a une fonction generatrice du bruit 
appartenant a un espace (Bj) prealablement defini de fonctions generatrices du bruit ; 

d) specification d'une norme ou d'une semi-norme dans Tespace des donnees ; 

15 e) specification d'une semi-norme dans I'espace des fonctions generatrices de bruit, 

pour laquelle chaque operateur de modelisation du bruit etablit sensibiement une relation 
isometrique entre T espace des fonctions generatrices du bruit et I'espace des donnees ; 

f) definition d'une fonction cout quantifiant I'ecart entre les mesures d'une part et Ja 
superposition de la reponse du modele et des bruits correles associes aux fonctions 

20 generatrices de bruit d'autre part ; et 

g) ajustement du modele et des fonctions generatrices du bruit en minimisant la 
fonction cout, par une methode algorithmique tirant paitie, des proprietes d'isometrie des 
operateurs de modelisation du bruit. 

2) Methode selon la revendication 1, dans laquelle on recherche la distribution en 
25 fonction de la profondeur de limpedance acoustique dans le milieu, les bruits correles 

affectant les donnees sont des ondes de tube identifiees chacune par des parametres 
caracterisant leur propagation, les donnees raesurees sont des donnees de PSV en utilisant 
des capteurs adaptes a detecter le deplacement des paiticules dans le milieu en reponse a 
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une excitation sismique localisee, remplacement des capteurs, la cluree d'enregistrement et 
les points d'echantillonnage en temps etant definis, et Toperateur de modelisation choisi 
associe les donnees synthetiques a une distribution d'impedance acoustique en fonction de 
la profondeur evaluee en temps de parcours et a la contrainte verticale mesuree en fonction 
du temps a la profondeur du premier capteur. 

3) Methode selon la reyendication 2, dans laquelle la fonction cout quantifiant 
Tecart est le carre de la semi-norme de cet ecait dans Tespace des donnees. 

4) Methode selon la revendication 2 bu 3, dans laquelle Tajustement du modele et 
des fonctions generatrices du bruit est obtenu par une methode de relaxation par bloc pour 
eliminer les inconnues conespondant a chaque fonction generatrice de bruits correles,\cette. 
methode de relaxation etant" mise en o^uvre a rinterieur des iterations d'un algorithme de 
quasi-Newton pour lecalcul du naodel^ . , . . . \ 

5) Methode selon Tune des revendications 2 a 4^ dans laquelle le calcul numerique^ 
de r image d'un modele par Toperateur de modelisation est effectue par resolution 
numerique de 1' equation des ondes ID pour )e rnodele considere en selectionnant des 
valeurs prises par le depl^ des: particules aux emplacements de capteurs ,et aux; 
points d'echantillonnage en templ^ prealablemerit specifies et en appliquant un.operateur 
propre a compenser les effets de divergence spheriqiie et d'attenuation. / ; 

6) Methode selon Tune des revendications 2 a 5, dans laquelle I'operateur 
numerique de, modelisation du bruit est un schdma numerique centre aux differences firiies 
pour discretiser reqiiaition de transport du bmit et la fonction generatrice de biuit 

\ intervenant conime Condition initiale au bord de la zone d'pbservdtion appartient a, un 
espace (Ej) de fonctions du temps a support dans un intervalle de temps donne. 

7) . Methode selon tiine des revendications 2 a 6, dans laquelle dans laquelle la 
semi-norme choisie pour I'espace des donnees est : 
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et la semi-norme choisie pour I'espace des fonctions generatrices de bruit est 


1 


\m\B = [AxAti (r-'' + r)' 


n=0 


8) Methode selon la revendication 1, dans laquelle on recherche la distribution de 
peiturbations, par rapport, a un modele de reference prealablement choisi, de I'impedance 

5 et de la vitesse dans la dite zone du milieu, les bruits correles affectant les donnees sont dus 
a des reflexions multiples dont la cinematique et les variations d' amplitude avec I'offset 
ont ete prealablement estimes, les donnees mesurees sont captees par des capteurs 
sismiques de surface, Templacement des dits capteurs, le mode d'excitation sismique, la 
duree d'enregistrement et le pas d'echantillonnage en temps etant definis, et Toperateur de 

10 modelisation etant defini via une linearisation de I'equation des ondes autour du modele de 
reference. 

9) Methode selon la revendication 8, dans laquelle la fonction cout quantifiant 
I'ecart est le carre de la semi-norme de cer ecart dans Tespace des donnees. 

10) Methode selon la revendication 8 ou 9, dans laquelle I'ajustement du modele et 
15 des fonctions generatrices du bruit est obtenu par une methode de relaxation par bloc pour 

eliminer les inconnues correspondant a chaque fonction generatrice de bruits correles, cette 
methode de relaxation etant mise en oeuvre a I'interieur des iterations d'un algorithme de 
gradient conjugue pour le calcul du modele. 
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S'il y a plus de trois inventeurs, utilisez plusieurs formulaires. Indiq'uez en haul a droite le N°.de la page suivi du nombre de pages. 


DATE ET SIGNATURE(S) 
DU (DES) DEMANDEUR(S) 
OU DU MANDATAIRE 
(Nom et qualite du signataire) 

Alfred ELMALEH, 
Chef du Departement Brevets 



La lol n°78-17 du 6 Janvier 1978 relative a rinformaticjue. aux fichiers et aux Itbertes s'applique aux reponses faites a ce formulaire. 
Elle garantit un droit d'acces et de rectification pour les donnees vous concernant aupres de TINPI, 


